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We use first principles density functional theory to calculate the electronic structure of the 
phenylthiolate (S-CeHs) self-assembled monolayer (SAM) on Cu(lll) and Au(lll) substrates. We 
find significant lateral dispersion of the SAM molecular states and discuss its implications for trans- 
port properties of the molecular wire array. We calculate the two photon photoemission spectra and 
the work function of the SAM on Cu(lll) and compare them with the available experimental data. 
Our results are used to discuss assignments of the observed spectral data and yield predictions for 
new electronic states due to the monolayer not yet accessed experimentally. 
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I. INTRODUCTION 

Recently there has been increased interest in the possi- 
bility of using organic molecules as electronic components 
in nanoscale devices. One of the most important issues is 
the role of molecule- metal contact. The charge injection 
at the metal-molecule interface is usually modelled by the 
energy dependent probability for a substrate electron to 
tunnel to the molecular wire. The magnitude of the tun- 
nelling probability depends on lineup of the molecular 
wire electronic states and the metal Fermi level. When 
a closely packed array of molecular wires is formed on 
the metal electrode, the molecular states form bands dis- 
persive in the lateral direction. This suggests that the 
tunnelling probability in addition to the energy of the 
injected electron depends also on the lateral component 
of the electron wavevector. 

The change of the metal work function due to a self- 
assembled Monolayer (SAM) may be employed in device 
applications to assist the charge injection into the molec- 
ular wire 1] . The work function of a crystal surface has 
generally two contributions ^]: (1) the electrostatic bar- 
rier due to the distortion in the charge distribution at 
the surface and (2) the many-body effects of the screened 
hole. It is important to know the relative contributions of 
the two components for the successful design of molecular 
electronic devices. 

In this paper we present Density Functional Theory 
(DFT) calculations of the electronic structure of SAMs 
of phenylthiolate (FT = -S-C6II5) on noble metal sub- 
strates Cu(lll) and Au(lll). We find significant lat- 
eral dispersion of the molecular electronic states derived 
from the molecular orbitals of the SAM. Our calcula- 
tions aim to address three essential questions: (1) what 
are the molecular states closest to metal Fermi level and 
therefore the most important for transport properties; 

(2) how large is the lateral dispersion of those states; 

(3) what is the largest contribution to the workfunction 
change associated with SAM formation, for which we find 



the electrostatic contribution from the polar molecules as 
determined by the dipole moment of the isolated array 
of molecules in addition to that of the bare metal sur- 
face. In addition we calculate two-photon photoemission 
(2PPE) spectra using DFT wavefunctions and the work 
function of the SAM/Cu(lll) system to compare with 
the experiments by Zhu et. al. 



II. METHOD OF CALCULATION 

We use the full potential linearized augmented plane 
wave (FP-LAPW) method with local orbital ex- 
tensions laj in the WIEN2k implementation 161]. The 
GGA |2i] exchange-correlation potential was used. Well- 
converged basis sets were employed with a 4.7 Ry plane 
wave cut off. Five special k-points were used to sample 
the two dimensional Brillouin zone. The surface plane 
was taken as the x, y plane. We used slab calculations 
based on a super cell with three metal layers covered on 
one side by a FT monolayer for electronic structure cal- 
culations, and six metal layers with both sides covered for 
work function calculations. The FT monolayer is of the 
so called \/3 x \/^R2iQP type [g], with one molecule per 
three substrate surface metal atoms. The former are ar- 
ranged in a triangular lattice with side \/3Rm-m , where 
Rm-m is the side of the triangular lattice formed by 
metal atoms M in the (111) plane. For M=Au and Cu, 
the corresponding cross-section area for molecule is, re- 
spectively, 21.6 A^ and 16.9 A^ HQ. The cross-sectional 
mean area of a FT molecule (when directed normal to the 
substrate) is 21.1 A^, slightly less than the above value 
based on the Au substrate, thus indicating the possibility 
of a small degree of tilt 8, 2\ , whereas no tilt is expected 
in the case of the Cu substrate. 

The geometry of the molecular array has been fully re- 
laxed on an unrelaxed substrate with metal atoms in ideal 
crystal positions. The monolayer on the Au substrate 
was found tilted by 18°, consistent with related experi- 



mental data |9| . In the calculation, the PT molecule was 
constrained so as to maintain C2V intramolecular symme- 
try jl3|- The lowest-energy binding site was the Au-Au 
bridge, analogous to the case of alkanethiols studied by 
Selloni et. al. Jj. The S atom height is 2.1 A. Starting 
with alignment normal to the substrate, with the projec- 
tion of the PT plane parallel to a nearest-neighbor Au-Au 
vector, the PT was tilted in the direction perpendicular 
to its plane (tilting in the orthogonal direction caused a 
sharp rise in energy). The calculated tilt relaxation en- 
ergy per molecule (0° to 18") was of the order of /c^T at 
room temperature, indicating the likelihood of apprecia- 
ble thermal fluctuations, similar to the situation noted in 

Si. 

In the case of the Cu substrate we find zero tilt angle 
for the monolayer due to the close spacing controlled by 
the Cu lattice constant, with the lowest energy when sul- 
fur binds to the hollow site. The height of the S atom is 
1.73 A. The binding energy for the bridge site is higher 
by 170 meV. The frozen lattice approximation used here 
may cause a shift of the energy difference in binding en- 
ergy of a few tens of meV for different adsorption sites, 
but it should not affect much the work function and the 
shape of the dispersion curves. 



III. BAND STRUCTURE RESULTS 

The electronic structure along the high symmetry di- 
rection of the 2D Brillouin zone is shown in Fig. |^) 
for the PT SAM on Cu(lll). In order to distinguish 
which bands originate from the SAM and which are due 
to the metallic substrate we calculate projected density 
of states of each band inside the spheres surrounding the 
atoms of the PT molecule 4, 5]. The Cu surface states 
were identified by the similarities of the band dispersion 
of the covered and clean Cu(lll) three layer substrates 
shown on Fig. Qd). At least seven molecular bands due 
to the SAM can be identified. These molecular bands 
can be characterized by (1) the energy at the F-point, (2) 
the effective mass for the F-M direction in cases where 
a parabolic fit of the dispersion in the vicinity of the F 
point is possible, (3) the bandwidth defined as the energy 
difference between the M and the F point, and (4) the 
major atomic contributions. These band characteristics 
are summarized in Tabled 

For the copper substrate the molecular levels closest 
to the Fermi level are the two HOMO's (bands 3 and 4), 
which are of S (3p) character. They have a very small 
dispersion (of about 0.2-0.3 eV) and the parabolic fits 
result in large effective masses (-2.0 me). The second oc- 
cupied state (SSo, at -0.9 eV) and the second unoccupied 
state (SS„, at 1.7 eV) are surface states, which arc at 
energies 0.0 eV and -1.6 eV, respectively, on the clean 
substrate in Fig. Qd). The discrepancy with the exper- 
imental value of the -0.4 eV surface state energy at F 
point is due to the finite size effect of our three layer slab 
geometry fl2|. The largest lateral dispersion is due to 
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TABLE I: Analysis of the electronic bandstructure for the 
PT SAM on Cu(lll), SAM on Au(lll), and the free molec- 
ular arrays at the same geometries. The energies of the 
bands are given relative to the Fermi level at the F-point 
{k — (0,0)). Effective masses are reported for the bands 
where a parabolic fit (along the F-M direction) is possible. 
The bandwidth is reported for bands where the energy dif- 
ference between the M and the F point can be identified (M 
denotes k = (Itt/Vs, 0)). The last column shows the band 
character for the SAM/substrate systems, where Cs stands 
for side phenyl carbons and Ct for next to sulfur and top 
carbons. 
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FIG. 1: (color online). Band structure results (energies rela- 
tive to Epermi). (a) PT SAM on 3 layers of Cu(lll). The size 
of the colored circles is proportional to the projected density 
of states of each band on a given atom of a molecule: sulfur 
- blue, Cs carbons - light blue, and Ct carbons - green, (b) 3 
layer slab of Cu(lll). From comparison of (a) and (b) four 
HOMO and three LUMO bands can be identified due to the 
SAM attachment. 
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FIG. 2: (color online). Band structure results (energies rela- 
tive to Epcrmi) (a) PT SAM on 3 layers of Au(lll). The size 
of the colored circles is proportional to the projected density 
of states of each band on a given atom of a molecule: sulfur 
- blue, Cs carbons - light blue, and Ct carbons - green, (b) 3 
layer slab of Au(lll). 
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FIG. 3: (color online), (a) Band structure of the free standing 
PT array, based on the Cu(lll) lattice constant and (b) a 
PT array tilted by 18°, based on the Au(lll) lattice constant 
(energies relative to Epcrmi). The size of the colored circles is 
proportional to the projected density of states of each band 
on a given atom of a molecule: sulfur - blue, Ca carbons - 
light blue, and Ct carbons - green. 



the bands 1 and 5 originating from the side phenyl car- 
bons (Cg for short). (The minimal separation between 
such atoms on neighboring PT molecules is 2.61 A.) The 
molecular bands 2, 6, and 7 are mainly due to the next 
to sulfur and top carbon atoms (C( for short). These Ct 
atoms form the longest intermolecular contacts and the 
lateral dispersion of the C(-type bands is not significant 
compared to that of side carbon Cg-type bands. There 
are certainly more electronic states in the band struc- 
ture due to the SAM, as is evident from the band count 
in Fig. (^, but these are high energy states (relative 
to Fermi level) and are not important for the conduc- 
tion properties of the monolayer. These states are also 
well delocalized in space such that projected density of 
states inside the muffin tin spheres does not give a unique 
identification of the molecular bands. As will be clear in 
Section Hvl the band 6 at 2.8 eV (the F point) and one at 
6 eV have the largest contribution to the 2PPE intensity 
and thus are candidates for the peaks observed by Zhu 
et. al. 0. 

When the PT is attached to the gold substrate, the 
bandwidths of all molecular levels are narrower than for 
the Cu case, because the nearest neighbor distance be- 
tween the molecules is larger due to the larger lattice 
constant of the underlying substrate. The Cs atoms on 
neighboring molecules of the SAM (18% tilted on Au) are 
separated by at least 3.13 A from each other, whereas 
Ct carbons are 5.0 A apart (i.e., V^Rau-Au)- The elec- 
tronic band structures of the SAM/Au(lll) and the clean 
Au(lll) surfaces are shown in Fig. Q. The band as- 
signments are summarized in table (QJ . As in the case of 
the SAM on Cu, the two nearly degenerate sulfur-based 
homo's (bands 3 and 4) are the molecular levels closest 
to the Fermi level. As the lattice constant increases (in 
proceeding from Cu to Au), the direct exchange mech- 



anism responsible for the negative effective mass decays 
faster than the substrate mediated mechanism (see be- 
low) . This causes the sign of the effective mass of one of 
the sulfur bands to become positive on the Au substrate. 
The most dispersive Cs bands 1 and 5 have bandwidth 
smaller by a factor of 2.5 than for the Cu substrate. The 
occupied surface state SSo is strongly hybridized with the 
sulfur 3p orbitals. 

The lateral dispersion has two contributions the direct 
overlap between the molecules and the substrate medi- 
ated interaction. To estimate the relative importance of 
the two mechanisms for the band dispersion we compare 
the lateral dispersions obtained for the PT arrays in the 
absence of the substrates (infinitely separated from the 
substrate) , but maintaining the same geometry as deter- 
mined with the substrate present. In this case the band 
structure shown in Fig. |21 has only a direct contribu- 
tion. The signs of the dispersions for bands 3 and 4 are 
negative and the bandwidth is smaller for the Au(lll) 
case (Fig. (Od)). In both cases the dispersion of the sul- 
fur bands is larger in the absence of the substrate. This 
suggests that substrate-mediated interactions have an ef- 
fect on the bandwidth opposite to that from the direct 
exchange. When the lattice constant is relatively large 
(i.e. the Au case), the substrate-mediated interaction 
even changes the sign of the band 4 dispersion. The cou- 
pling of the Cs atoms to the substrate is much smaller, 
so that substrate-mediated interaction in bands 1 and 5 
is responsible for only 10% of the total bandwidth. 



IV. 2PPE SPECTRA 

The electronic structure of the organic SAM can be 
directly probed in the 2PPE experiment M^. There are 
three electronic states involved in the 2PPE spectra: ini- 



tial (z), intermediate (k), and the final (/). If one of those 
states is locahzed and has the largest dipole matrix el- 
ement in comparison with the rest of the states, then 
three scenarios are possible for the variation of kinetic en- 
ergy {7Ekm) with respect to the excitation photon energy 
hojL , which is chosen to be less than the workfunction of 
the substrate: (1) lEkin ~ '2fiLUL; (2) ?-E/cm ~ ^l, and 
(3) invariance with respect to Hlul. For PT on Cu(lll), 
two cases are reported Q]: case (2) in which the inter- 
mediate state (i) is assigned as a LUMO lying below the 
vacuum at 3.3 eV and case(3), in which the final state (f) 
is assigned as a LUMO lying above the vacuum at 6.4 eV 
(relative to Epormi)- Case (1), corresponding to an initial 
HOMO state and was not observed in Q- 

In the present work we use DFT wavefunctions to cal- 
culate the 2PPE spectra at the P-point. We find two 
states to contribute most to the 2PPE intensity: a LUMO 
(band 6) at 2.8 eV (of Ct carbon character) and high en- 
ergy bands at around 6.0 eV, whose character is complex 
due to hybridization of molecular and Cu states. This 
assignment is consistent with the conclusion drawn in [3| 
based on similarities of spectra taken for SAMs based on 
PT and alkanethiolates of different lengths. 

There are two possible 2PPE excitation mechanisms: 
direct and indirect. In the direct process with one-color 
pump and probe laser frequencies, the 2PPE intensity for 
a- light polarization is [iJi UM • 

2 

< f\pa\k >< k\pa\i > 



E 

i<k<f 



Efc - Ei - hiUL + iPfe 



Mi-fk) (1) 



where Pfc is the lifetime broadening of intermediate level 
k and a denotes the Cartesian component (x,y, or z). The 
sum over intermediate states in Eq. Q includes inter- 
ference effects. The delta function may be replaced by 
a Lorentzian to mimic the lifetime of the final state and 
the width of the laser pulse. The result for the P-point is 
shown in Fig. 2|where intensity contributions Eq. (Q are 
summed over all initial states i. The calculated curves 
were convoluted with a Lorentzian function (P = 0.1 eV) 
to account for the spectrometer resolution. We use in- 
termediate and final level lifetimes P/c — Tf = 0.1 eV. 
There is one pronounced peak at about 6.0 eV, whose 
intensity is resonant with the laser frequency at about 
3.6 eV. As the excitation frequency increases the inten- 
sity transfers to higher energy final states (~ 7 eV) . The 
resonance enhancement factor is sensitive to the choice 
of the damping parameters, and even for the large value 
(0.1 eV) used in Fig. Q, the intensity of the out-of-plane 
polarization varies by an order of magnitude. The reso- 
nance effect is related to the large dipole matrix element 
between the occupied surface state SSo and the second 
LUMO (band 6) separated by 3.7 eV. For the in-plane po- 
larization there is no strong laser frequency dependence 
of the 2PPE intensity, because there is no single pair of 
states which dominates the sum in Eq. Q. Therefore 
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FIG. 4: Prediction of the 2PPE spectra for the direct mech- 
anism according to Eq. (0 for (a) out of plane and (b) in 
plane light polarizations. 



the ratio of the out-of-plane and the in-plane polariza- 
tions in the direct mechanism model should be resonant 
with the laser frequency. 

In the indirect mechanism the intermediate level can 
be populated by some incoherent process, for example, 
by the tunnelling of the photoexcited electron from the 
substrate to the molecule. In the second step, a photon 
promotes the photoexcited electron to the vacuum with 
probability given by the one photon absorption crossec- 
tion: 



l'^''"{k, /, lol) - SihiOL - Ey + Efe) |< f\p^\k >f 

(l-/fc)(l-//)/fc(EFcrmi + fi^L) (2) 



where Fermi-Dirac factors / require that the intermedi- 
ate and final states are empty and the intermediate state 
is less than one quantum of photon energy above the 
Fermi level EFormi- It is assumed that all the intermedi- 
ate states k have the equal lifetimes and probabilities to 
be excited in the first step. The result for the two polar- 
izations and the same damping parameters as in direct 
process is shown in Fig. [3 There are no pronounced res- 
onances and the ratio of intensities for the out-of-plane 
and the in-plane polarizations is relatively independent of 
the laser frequency. Namely, the ratios of the integrated 
peak intensities at 6 eV for ujl =3.4, 3.5, 3.7, and 3.9 are 
1.8, 4.6, 3.8, and 4.0 respectively. This is consistent with 
the experimentally reported value of w 4 at lul — 3.7 
eV ^. Comparison of our analysis with the available 
experimental data 3] leads us to suggest that the most 
likely mechanism for the two photon photoemission pro- 
cess is indirect. An enhancement of 2PPE intensity with 
the laser frequency was observed for SCH3 on Ag(lll) 
by the Harris group [Ig . This is may be the same effect 
we find here for the direct mechanism. 
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FIG. 5: Substrate assisted 2PPE spectra ( indirect mecha- 
nism) calculated with Eq. ||5J for (a) out of plane and (b) in 
plane light polarizations. 
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FIG. 6; (color online). Coulomb potential averaged in the 
lateral direction versus z-direction for (a) the clean Cu(lll), 
(b) the single monolayer, and (c) the SAM on Cu(lll). The 
dipole moments according to Eq. l|HJl are shown by the arrows. 



V. WORK FUNCTION CALCULATION 

For the work function calculations we use a six Cu 
layer substrate covered with a SAM on both sides. The 
inversion symmetry of the supercell with two equivalent 
surfaces forces the net dipole and the electrostatic poten- 
tial drop in the vacuum region to be zero. To separate 
the electrostatic and exchange correlation contributions 
to the work function we calculate the averaged surface 
dipole moment by integrating over the net charge den- 
sity of the half unit cell: 



P 



Z/2 



z[p(a;, y, z) - p+{x, y, z)]dxdydz, 



(3) 



where p and p+ are valence electron and ion charge den- 
sities respectively. The planes z = and z = Z/2 contain 
the inversion centers of the supecell. The Coulomb po- 
tential averaged in the lateral direction is shown on Fig. |B1 
along with that for the clean substrate and the two SAM 
arrays in the absence of the substrate for the same geom- 
etry. The work function of the clean Cu(lll) is found to 
be 5.0 eV, in good agreement with the experimental value 
of 4.9 eV. The work function of the covered substrate is 
reduced to 2.7 eV compared with the experimental value 
of3.7eV 0. 

The isolated SAM layer has surface dipole Psam — 
1.19 (Debye/molecule) , and the expected Coulomb po- 
tential drop of 47rPsAMM = 2.67 eV (where A = IG.gA^ 
is the lateral area per molecule) agrees well with the ac- 
tual result 5V = 6.5-3.8 = 2.7 eV (see Fig. (0d)). Anal- 
ogously we define the electrostatic contribution to the 
clean Cu surface workfunction to be AttPcu/A = 4.6 eV 
(see Fig. (|2K))- This leaves a 0.4 eV contribution due 
to exchange-correlation effects. In the case of the SAM 
on Cu, the net surface dipole is approximately equal to 
the sum of dipoles of the clean Cu surface and the iso- 
lated SAM array -2.06 -I- 1.19 = -0.87 Debye. The 



difference between the simple estimate and the actual 
result (—0.076 Debye) is due to the charge transfer at 
the interface. The laterally averaged charge densities are 
shown in Fig. (|7^). The charge transfer electron den- 
sity shown in Fig. ^7}p) is very small (less then 0.1 e 
per molecule). We calculate the workfunction change of 
the covered substrate to be -2.3 eV, where -2.8 eV is the 
electrostatic and 0.5 eV is the exchange-correlation con- 
tribution, respectively. 

As we have shown the electrostatic contribution to the 
workfunction change can be immediately estimated from 
the dipole moment of the isolated SAM array. However, 
the isolated SAM dipole moment itself is not equal to 
that of isolated molecule. We have used the molecular 
NRLMOL code jlS:] and find the dipole moment of the 
S-CgHg molecule to be 3.33 Debye. This value is three- 
fold larger than the dipole per molecule in the array, 1.19 
Debye. We have also performed constrained calculations 
by forcing occupations of two nearly degenerate HOMO 
orbitals: Px (in the molecular plane (xz)) and Py (perpen- 
dicular to the molecular plane). The unconstrained spin 
restricted ( ie, pure doublet) calculation results in 1.88 
electrons in px and 1.12 electrons in py orbitals, with a 
net dipole moment of 3.33 Debye. By forcing the Py or- 
bital to be fully occupied the net dipole moment drops to 
1.24 Debye, very close to the SAM value. This suggests 
that the charge redistribution between the two HOMO 
orbitals in the SAM is the primary effect on the net dipole 
moment. 



VI. DISCUSSION AND CONCLUSIONS 

We have calculated the electronic structure of the 
phenylthiolate (PT=S-C6H5) SAM on Cu(lll) and 
Au(lll) substrates. We identified four HOMO's and 
three LUMO's for each system. A sulfur- type HOMO 
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FIG. 7: (color online), (a) The laterally average charges den- 
sities versus z-direction for the clean Cu(lll) (blue) - pi, 
isolated SAM array (red) - p2 and two-sided SAM/Cu(lll) 
(black) - p3. (b) The laterally averaged charge density differ- 
ence P3 — pi — p2 versus z-direction which integrates to the 
induced dipole P=0.08 (D) due to the charge transfer. 



(band 4) is the closest to the Fermi level for both sub- 
strates. It has two competing contributions to the lateral 
dispersion: (1) the direct exchange between the neighbor- 
ing orbitals leads to a negative effective mass, whereas (2) 
the substrate-mediated interaction gives a positive effec- 
tive mass contribution. The overall sign of the effective 
mass is changed from a negative to a positive value in go- 
ing from the Cu to the Au substrate. The most dispersive 
molecular bands are of 'side carbon' character (Cg). The 
significant lateral dispersion of the electronic states de- 
rived from molecular orbitals of the SAM suggests that 
the charge injection probability depends on the lateral 
component of the wavevector as well as the energy of the 
injected electron. This implies that the conductance of 
N molecules in an array is not simply N times the con- 
ductance of one, but further theoretical and experimental 
work is needed to quantify the lateral dispersion effect on 
transport. 

We have calculated the 2PPE spectra of a PT SAM 
on Cu(lll) using the DFT wavefunctions. In the di- 
rect mechanism we predict a very strong resonance of 
the 2PPE intensity for out-of-plane polarization. For the 
in-plane polarization only small intensity variation with 
laser frequency is found. For the indirect mechanism we 
find a moderate intensity dependence on the excitation 
photon energy for both polarizations with a nearly con- 
stant ratio of the out-of-plane to the in-plane signals of 
about 3-4 , consistent with related experimental rsults 



Q. This fact suggests that the indirect mechanism con- 
tributes the most intensity to the 2PPE resonance ob- 
served in Q ( designated as resonance A). The energy 
positions of the observed molecular states at 6.4 eV and 
3.4 eV are in reasonable correspondence with our cal- 
culated values of 6.0 eV and 2.8 eV, respectively. These 
molecular states are mainly due to LUMOs dominated 
by the PT carbon atoms adjacent to the S atoms. 

Calculation of the workfunction for the PT-SAM on 
Cu(lll) yields a value of 2.7 eV, which may be compared 
with the experimental value of 3.7 eV. To reconcile the- 
ory and experiment one may assume a lower coverage of 
the SAM in the experimental system. The larger effec- 
tive area per molecule reduces the surface dipole moment 
and hence the Coulomb contribution to the workfunction 
change. According to our calculations one has to assume 
A K, 23.8 A^ to agree with experiment. This estimate is 
much larger than that for the \/3 x \/3^30° structure, 
which implies A = 16.9 A^. Alternatively, the interpre- 
tation of the workfunction measurements in one photon 
photoemission spectra can be complicated if there is an 
exit barrier for substrate electrons to tunnel through the 
interface. Then the onset of the one-photon photoemis- 
sion signal may be due to the finite energy electrons with 
respect to the vacuum level. This in turn would imply 
that the agreement between the theory and the exper- 
iment for the energies of the LUMO orbitals measured 
in 2PPE is accidental. Indeed DFT is known to under- 
estimate the bandgaps in semiconductors. On the other 
hand the energy differences between the empty states are 
more reliable in DFT and in the experiment the energy 
difference between the two LUMO's is also independent 
of the workfunction value. We find an electrostatic con- 
tribution to the work function change of 2.8 eV ( asso- 
ciated with SAM formation) to be the dominant effect. 
Therefore the dipole moment of the molecular SAM in 
the absence of the metal substrate can be used to predict 
the the workfunction change of the covered substrate pro- 
vided the charge transfer and correlation effects do not 
contribute much to the workfunction change as in the 
case studied here. 
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